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IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 

In re Patent Application of 

CHAPMAN, Barry L.W. Atty. Ref.: 604-632 

Serial No. Unassigned Group: 
Filed: March 29, 2002 Examiner: 
For: PERMANENT MAGNET AND SHIM DESIGN 

*********** 

March 29, 2002 

Assistant Commissioner for Patents 
Washington, DC 20231 

Sir: 

PRELIMINARY AMENDMENT 

Please amend the above application as follows: 

IN THE CLAIMS 

Please substitute the following amended claims for corresponding claims 
previously presented. A copy of the amended claims showing current revisions 
is attached. 

6. (Amended) A method according to claim 1 wherein the desired field is 
produced by a distribution of magnetic material over one or more flat surfaces. 

615806 



CHAPMAN, Barry L.W. 
Serial No. Unassigned 

7. (Amended) A method according to claim 1 wherein the desired field is 
produced by a distribution of magnetic material over an elliptical cylindrical shell. 

8. (Amended) A method according to claim 1 wherein the desired field is 
produced by a distribution of magnetic material over the surface of a circular 
cylindrical shell. 

9. (Amended) A method according to claim 7 wherein the thickness of 
magnetic material is varied in an azimuthal direction in accordance with the 
desired field. 

12. (Amended) A method according to claim 1 wherein an additional, 
hypothetical, arbitrary distribution of magnetic material over said surface is 
incorporated into said set. 

13. (Amended) A method according to claim 2 comprising: 
-calculating the distribution of magnetic material to produce an optimum 

approximation to a desired magnetic field in the ROT, in accordance with the 
method described herein 

-determining the total resultant field in the ROI taking into account the 
thicknesses and placements of magnetic materials involved 
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-calculating a corrected magnetic material distribution representative of the 
difference between said total resultant field and said desired field 

-repeating tlie latter two steps until the field attained approximates the 
desired one in the ROT to a sufficient degree. 

1 4. (Amended) A method according to claim 1 where the magnetic fields of 
said set of harmonics are evaluated by an actual step involving numerical 
integration in real space. 

15. (Amended) A method according to claim 1 where the Fourier 
transformers of the magnetic fields of said set of harmonics are evaluated directly 
in Fourier Kernel space. 

16. (Amended) A permanent magnetic structure designed in accordance 
with the method of claim 1 . 
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REMARKS 

The above amendments have been made to place the application in a 
more traditional format. Attached hereto is a marked-up version of the changes 
made to the claims by the current amendment. The attached pages are 
captioned " Version With Markings To Show Changes Made " 

Respectfully submitted, 
NIXON & VANDERHYE P.C. 

Rec/. No. 29,009 

11 00 North Glebe Road, 8th Floor 
Arlington, VA 22201-4714 
Telephone: (703) 816-4000 
Facsimile: (703)816-4100 



By: 



LCM:lks 
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VERSION WITH MARKINGS TO SHOW CHANGES MADE 
IN THE CLAIMS 



hi 



m 



mi 



6. (Amended) A method according to [any of claims 1 to 5] claim 1 
wherein the desired field is produced by a distribution of magnetic material over 
one or more flat surfaces. 



^ 7. (Amended) A method according to [any of claims 1 to 5] claim 1 

wherein the desired field is produced by a distribution of magnetic material over 
an elliptical cylindrical shell. 

Ill 



8. (Amended) A method according to [any of claims 1 to 5] claim 1 
9 wherein the desired field is produced by a distribution of magnetic material over 
the surface of a circular cylindrical shell. 



9. (Amended) A method according to claim 7 [or claim 8] wherein the 
thickness of magnetic material is varied in an azimuthal direction in accordance 
with the desired field. 



12. (Amended) A method according to [any of claims 1-1 11 claim 1 
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Wherein an additional, hypothetical, arbitrary distribution of magnetic material 
over said surface is Incorporated into said set. 

13. (Amended) A method according to [claims 2 or 3] claim_2 comprising: 
-calculating the distribution of magnetic material to produce an optimum 

approximation to a desired magnetic field in the ROT, in accordance with the 
method described herein 

-determining the total resultant field in the ROI taking into account the 
thicknesses and placements of magnetic materials involved 

-calculating a corrected magnetic material distribution representative of the 
difference between said total resultant field and said desired field 

-repeating the latter two steps until the field attained approximates the 
desired one in the ROT to a sufficient degree. 

14. (Amended) A method according to [any of claims 1 to 5] claim 1 where 
the magnetic fields of said set of harmonics are evaluated by an actual step 
involving numerical integration in real space. 

15. (Amended) A method according to [any of claims 1 to 5] claim 1 where 
the Fourier transformers of the magnetic fields of said set of harmonics are 
evaluated directly in Fourier Kernel space. 
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16. (Amended) A permanent magnetic structure designed in accordance 
with tlie method of [any of claims 1 to 15] claim 1. 
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PERMANENT MAGNET AND SfflM DESIGN 

FIELD OF THE INVENTION 
This invention relates to the design of permanent magnetic stractures for 
5 generating a specific magnetic field in a Region Of Interest (ROI) or modifying an 
existing magnetic field in an ROI using shims* 

BACKGROUND OF THE INVENTION 
Permanent magnets are of great value in many fields of science and 
engineering. In particular, they find application in the field of Magnetic Resonance 
10 Imaging (MRI) and Spectroscopy (MRS). Within these latter apphcations it is of 

paramount importance that extremely high degrees of magnetic field homogeneity 
(of the order of one part per million) are achieved within an ROI. 
In permanent magnets the distribution of the magnetic material determines 
how the resulting magnetic field varies over space. In both permanent magnetic and 
15 electromagnetic MRI and MRS machines, manufacturing imperfections result in 
unacceptable inhomogeneities. To improve the homogeneity of an existing field 
over an ROI it is known to place permanent magnetic materials (termed passive 
shims) in the magnetic field to achieve this improvement of homogeneity. 

The problem of distributing magnetic material to produce, or modify an 
20 existing magnetic field is one of a class known collectively as the near field 
problem. In particular, it requires the determination of the distribution of material 
to produce a field of particular character - varying in a specific manner in space, for 
example a imiform field in a proximate region. The total magnetic field fi^om a 
magnetic structure is a nonlinear convolution of the spatial distribution of material 
25 with the field of a magnetic dipole. While it is relatively straightforward to 
determine the resulting field fi:'om a given distribution of magnetic material, 
generally, the components of a convolution are not invertable. Consequently, 
detemiining the optimum distribution of magnetic material to produce a magnetic 
field of a specific type in an ROI has proven to be extremely difficult. 
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One known method of electromagnetic coil design is the Target field approach 
as disclosed in European Patent No. EP-A- 0 252 634 or its US equivalent No. 
4896129. As it is based on an assumed a priori mathematical expression for the 
total field, an optimum resulting field is hard to achieve. 
5 An alternative approach is disclosed in European Patent EP-A-0601 101 and its 

US equivalent No. 5266913, the content of which is hereby incorporated by 
reference. It is concerned with the design of shielded electromagnetic coils, and 
employs forward substitution and least squares minimization to optimize the field 
produced by electromagnets specifically within the ROI, whilst achieving a null 
10 field outside of the confines of the coil structure. 

There are significant differences in the form of the magnetic fields generated 
by active electromagnets, and those resulting fi-om the use of magnetic materials. 
Consequently, electromagnetic coil design and the design of permanent magnets 
involve quite separate design considerations. 
15 The object of the present invention is to achieve an improved method of 

designing permanent magnetic structures. 

SUMMARY OF THE INVENTION 
In accordance with the present invention there is provided a method of 
designing a permanent magnetic structure for generating a desired magnetic field in a 
20 region of interest (ROI) comprising: 

-setting the dimensions of the permanent magnetic structure 
-setting a hypothetical distribution of magnetic material over the surface of 
said structure for each member of a finite set of spatial fi-equency harmonics of unit 
magnetisation amplitude 
25 -calculating the magnetic field in tfie ROI ft^om each member of said set 

-calculating amplitudes of said harmonics by the method of least squares in 
order to achieve the optimum approximation to the desired field in the ROI 

-scaling the members of said set according to said optimised harmonic 
amplitudes 



-2- 



wo 01/29570 



PCT/GBOO/03905 



-establishing, by summing said scaled members, the required distribution of 
magnetic material in order to generate the desired magnetic field. 

The permanent magnetic structure may either be a pemianent magnet 
assembly for producing a magnetic field ab initio, or may be a permanent magnetic 
5 stmcture for modifjdng an existing magnetic field, by means of shimming. 

The method may also be applied iteratively in order to take into account the 
finite thickness of the required distribution magnetic material, which will inevitably 
have influence on the field. 

Further preferred embodiments of the present invention are presented in the 
1 0 dependent claims. 

BRffiF DESCRIPTION OF THE DIAGRAMS 
Figure 1 shows a schematic representation two parallel, finite, flat plates upon 
which magnetic material is distributed. 

Figure 2 shows an alternative distribution composed of a hollow elliptical, 
15 cylindrical stracture upon which magnetic material is distributed. 

Figure 3 shows the magnetic dipolar field (Bd) produced by magnetic material, 
resolved into orthogonal components Bdr , Bdz and Bde. 

Figure 4 shows schematically a hollow, cylindrical structure with a circular cross- 
section, upon which magnetic material is distributed. 
20 Figure 5 shows the pemiissible, even, cosinusoidal axial harmonics of the magnetic 
material distribution on a structure of finite length. 

Figure 6 is a plot of the axial magnetic field Brz against distance z resulting fi-om an 
infinitesimal circular ring of imiformly magnetic material of unit magnetization 
calculated directly by numerical integration. 
25 Figure 7 shows a plot of the Fourier Transform (FT) in the axial direction (z) of the 
axial magnetic field of Figure 6. 

Figure 8 shows a plot of the FT, in the axial direction, of the axial magnetic field 
resulting firom an infinitesimal circular ring of uniformly magnetic material of imit 
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magnetization evaluated by an alternative, analytical method, directly in the Fourier 
Kernel domain (k-space). 

Figure 9 shows a plot of the axial distribution of magnetic material that will achieve 
an optimal, uniform field, designed in accordance with a first embodiment of the 

5 present invention. 

Figure 10 shows a schematic representation of a magnetic structure designed in 
accordance with the invention, to generate the distribution shown in Figure 9. 
Figure 1 1 shows a plot of the resulting axial magnetic field in the axial direction 
produced by the distribution of magnetic material shown in Figure 9. 

10 Figure 12 shows a plot of the deviation of the above magnetic field fi-om the 
desired, uniform one, in parts per milhon (ppm), over the ROI. 
Figure 13 shows a plot of the deviation of the above field over, and proximate to, 
the ROI. 

Figure 14 shows how a distribution of paramagnetic and diamagnetic materials on a 
15 surface can be arranged to produce a sinusoidal variation of magnetization along 
that surface. 

Figure 15 shows a plot of the amplitude of the optimum axial distribution of 
magnetic material over a cylindrical structure to provide a magnetic field that varies 
linearly in a transverse direction in the ROI, designed in accordance with the present 
20 invention. 

Figure 16 shows a schematic representation of how a shim magnetic structure 
designed in accordance with the invention and using diamagnetic and paramagnetic 
materials as illustrated in Figure 14 could be employed to generate the distribution 
shown in Figure 15. 

25 Figure 17 shows a plot of the resulting axial component of the gradient magnetic 
field at a series of radii in the axial direction produced by the distribution of 
magnetic material shown in Figure 15. 
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Figure 18 shows how tiie sinusoidal variation of thickness in a single unidirectional 
magnetic material can provide a sinusoidal variation of magnetization for the 
purpose of shimming. 
Theory. 

5 The design of magnets for MRI and MRS is commonly based upon the forms 

of simple geometric surfaces. The two most common of these are: a stmcture 
composed of parallel, finite, flat surfaces between which a uniform magnetic field is 
created as illustrated in Figure 1; and a structure with an elliptical, cylindrical 
opening within which a uniform field is produced by the surrounding distribution of 

1 0 material as illustrated in Figure 2. 

With reference to Figure 3, the magnetic dipolar vector field (Bd) produced by 
magnetic material can be resolved into its cylindrical polar components Bdz^Bde and 
Bdr. In these coordinates (r, 0, z), the vector field of such an infinitesimal dipole 
has the general form: 

15 

Bd,(r,z) = m (2z^ •?)IA7t (r^ + z^^ (la) 
Bd0{r,z) = 0 (lb) 

and 

20 Bdr(r,z) = ^ m 3zr / 471 (r^-h z^^ (Ic) 

where ki is the permeability of the medium in the ROI and m is the magnetic dipole 
moment. The field due to the circulation of an electron in an atom ^proximates 
this very closely for distances (r^ + z^)^^^ that are very much greater that the radius 
25 of the orbit (typically 10"^^ m). For the purposes of passive magnet design within 
the present context this approximation holds extremely well. 

The magnetic field firom matter in bulk can be derived from equation 1. For 
the axial component of the magnetic field this is given by: 
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Bz(r,z) = ii J {M,(x,y,z) (2(z-2f - [(x^xf + (y-yf ]) / 47c|r - rf } dx' (2) 

where the underlimng of individual items indicates that they are vector quantities 
and Mz is the axial component of magnetization (M). Similar expressions can be 
5 derived for the field in the other directions. 

For paramagnetic and diamagnetic materials the magnetization M is given by: 

M = X„,Bm (3) 

10 where Bm is the main, magnetizing field and Xm is the susceptibility. In general, 
Xm is a tensor of rank two. However, for isotropic materials it reduces to one of 
rank zero, e.g., for diamagnetic copper Xm = -LO x 10'^; and for paramagnetic 
aluminium Xm = 23 x 10'^. 

Equation (2) shows that the magnetic field is the convolution of the 

15 magnetization distribution with a spatial dispersion function. Problems of this 
nature can be conveniently solved by Fourier transforming into the kernel space (k- 
space) where the convolution becomes a simple multiplication, which greatly 
facilitates subsequent calculations. Indeed, for certain simple geometric structures, 
such as cylinders, there is no simple analytic expression for the magnetic field in 

20 real space, while there is one for the field in k-space. The existence of such an 
expression in k-space expedites the process of handling the above equations but it is 
not a necessary condition of the design method of the present invention. Indeed, the 
method is generally applicable, as is illustrated in the following examples. 

25 DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The design of a permanent magnetic structure according to the invention will 
now be described in general terms* The geometry and overall dimensions of the 
desired structure are set, as is a hypothetical distribution of magnetic material over 
the surface of the structure for each member of a finite set of spatial fi^equency 
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harmonics of unit magnetization amplitude. The field within the ROI from an 
element of the magnetic material of unit magnetization can then be determined. 
When mmierical methods are employed the field must be determined to a point 
where it approximates zero in the directions in which it is to be Fourier Transformed 

5 (FT'ed). This field is then FT'ed. The resulting field in Fourier kernel space (k- 
space) is multipUed by the FT*s of a set of the permissible harmonics (the lowest 
spatial frequency harmonics of the initial magnetic field distribution) of the 
magnetization on the surface. These resulting fields in k-space are subsequently 
inverse FPed to yield the fields in real space that are produced by each of the 

10 permissible harmonics in the selected set. These fields are fitted to the required 
field within the ROI to determine the amplitudes of the magnetization harmonics 
that minimize the sum of the squares of the deviation of the actual field from the 
required one. Fitting is performed using sets containing successively higher order 
harmonic components until the sum of the squares of the deviations of the 

15 achievable field from the desired one becomes sufficiently small. Eventually a 
point will be reached where deviations due to manufacturing tolerances will 
dominate, and no fiarther practical improvement will be possible. 

For simple geometric structures, such as a circular cylinder, this process can 
be facilitated by directly calculating the fields in Fourier k-space. 

20 The advantage of least squares fitting lies in its convenient mathematical 

properties. For example, it commutes with the Fourier transform; as the FT itself is 
a specific example of least squares minimization* 

The following discussion serves to illustrate the method of the present 
invention only and is not intraded to Kmit its gmeral applicabiUty. 

25 In a firat embodiment of the present invention, the method is applied to design 

a magnet where magnetic material is distributed on the curved surface of a cylinder 
to produce a uniform field. 

For the present purposes it will be supposed that it is required to create a 
uniform axial magnetic field Bz of 0.1 T over a cylindrical ROI of length 0.4m and 
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radius 0.125in by distributing magnetic material on the curved surface of a circular 
cylinder, such as that shown in Figure 4, having a length L of 2m and radius a of 
0.25m. The cylindrical ROI and surrounding cylindrical structure are concentric 
and coaxial along the z axis. For simplicity it will be assumed that the thickness of 

5 the magnetic shell is small compared to the diameter of the magnet, for the present 
purposes infinitesimal. In practice, a more elaborate calculation may be required to 
take into account the thickness, as described later. 

Conveniently, optimization is performed over a representative set of radii 
spanning the ROI, e.g. six radii at intervals of 20% from 0 to 100% of the radius of 

10 the ROI within the ROI However, for the present example, optimization over just 
two curved cylindrical surfaces of diameter of Om and 0.25 m and lengths 0.4m, was 
sufficient, as demonstrated by table 3 and Figure 12 (explained below). 

The distribution of magnetic material on the surface of the structure shown 
in figure 4 can be conveniently described in terms of a seri^ of orthogonal 

15 azimuthal and axial, sinusoidal and cosinusoidal, harmonics. The azimuthal 
harmonics form a discrete set of cosinusoidal terms cos (m(|> + a) and the axial 
harmonics form a truncated set of H(z) cos (nO + p), as illustrated in Figure 5, 
where H(z) is the unit step function given by: 

20 forlzl <L/2 H(z)=l 

(4) 

forlzl>U2 H(z) = 0 

where L is the length of the cylindrical surface and a and p are phase terms that take 
25 account of both the sinusoidal and cosinusoidal tmns. 

In many applications the set of harmonics is severely curtailed by 

considerations of symmetry. For example, for the present uniform field there is no 

azimuthal variation and all azimuthal terms other than m=0 are zero. 

Consequently, the problem of optimization reduces to a one-dimensional one of 
30 determining just the axial variation. Further, as the ROI and the structure are 
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concentric the resulting even symmetry of the system implies that all sinusoidal 
terms are also zero leaving only the purely cosinusoidal terms. It remains to 
detennine the amplitudes of each of the significant cosinusoidal terms to produce an 
optimum uniformity of the field within the ROI on the basis of the minimum sum of 

5 least squares deviations of the field. 

Now the field firom an elemental magnetic dipolar ring of unit magnetization 
can be determined by numerically integrating equation (2). Figure 6 shows how the 
resulting magnetic field varies with distance, in the axial direction, at a radius of 0.1 
m from such a dipolar ring of diameter 0.5 m. 

10 From equation (2) we also know that the magnetic field fi*om individual axial 

harmonics is the convolution of the field firom this ring with the harmonic. This 
convolution is most conveniently performed by simple multiplication in the Fourier 
domain. Consequently, the next step in the design process is to FT this field in the 
axial direction z; the resulting FT is shown in Figure 7. 

15 Alternatively, for a simple geometric surface such as this, the FT in the axial 

direction of the axial component Bxz of the ring field, Br, can be expressed 
analytically as: 

for r < a Br^ (r,m,k) - - n a Ulkj r) Km(|k|a) / 2n (5a) 
20 for r > a Brz (r,m,k) = - jx a k^ In.(|k| a) Kn,(lkir) / 2n (5b) 



where Im and are the modified Bessel fimctions. 

Figure 8 shows the FT in the axial direction of the axial component Brz of the 
ring field, Br, evaluated directly by equation 5a. This is substantially identical to the 
25 field derived by numerically integrating the generally applicable equation (2) and 
then Fourier Transforming it, as can be seen by comparing Figures 7 and 8. 

Although evaluating the FT directly by means of equation (5a) facilitates the 
process, this mefliod is not applicable to structures where an exact solution in the 
Fourier kernel domain is not available. The longer route of numerical integrating 
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equation (2) and subsequent FT'ing has general application for any structure, as it 
does not rely upon there being an analytic form of the equations in either real or k- 
space. 

Now, for an open cylinder of restricted length, L, the total magnetization, Mz, 
5 can be expressed as the sum of the permissible axial hamionic modes, n, on the 
surface as follows: 

Mz(a,z) = S„CnH(z)cos(2 7cnz/L-<j>) (6) 

10 where Cn is the amplitude of the nth magnetization harmonic. The FT of diis in the 
z-direction is: 

Mz(a, k) - 7c L Sn Cn [e^* sinc(k L/2 - nn) + e"^^ sinc(k U2 + nn)] (7) 

15 

Therefore, the FT of the total resultant field from the cylindrical surface is 
simply the product of the FT of the magnetization with the FT of the elementary 
ring field within this domain, i.e.: 

20 Bz(r, m, k) = M^Ca, k) Brz (r, m, k) (8) 

Optimization of the distribution can is then solved by least squares 
minimization of the first few terms (Cn) over the ROI in real, as against kernel, 
space. The number of terms is dependent upon the geometries of the coil structure 
25 and the ROI This can be determined by solving with increasing numbers of terms 
(Cn) until a sufficiently homogeneous field has been attained. 

In the above example the first 3, 4 and 5 harmonic amplitudes (coefficients) 
w^e employed and the worst deviation firom the field in the ROI determined (see 
Table 1). Optimization over more than 5 harmonic coefficients has no further 



wo 01/29570 



PCT/GBOO/03905 



significant effect. As can be seen from Table 1, the maximum deviation of the field 
in the ROI drops to - Ippm within the first 5 terms. In any case, practical 
considerations, such as manufacturing tolerances and the finite thickness of the 
structure, mean that fiirther improvements in field uniformity are impractical. 
5 Further approximation can be achieved by shimming the magnet in some fashion, 
for example along the lines discussed in the second embodiment below. 



Table 1 : Field deviation with number of harmonic components employed. 



No. of Harmonics 


Range of Harmonics 


Max deviation in ROI 


3 


0-2 


200 ppm 


4 


0-3 


16ppin 


5 


0-4 


1 ppm 



10 The harmonic amplitudes of table 2 were used to determine the axial 

distribution of magnetic material (M - 1.28 x 10^ amp/m, equivalent to the 
magnetization of saturated iron) on the surface of the hollow cylindrical structure, as 
shown in Figure 9. 

Figure 10 shows a practical implementation of a magnetic stmcture designed 
1 5 in accordance with the invention, to generate the distribution shown in Figure 9. 



Table 2: The harmonic amplitudes for a uniform field. 



Harmonic No. 


Amplitude (amps/m) 


0 


-1,143,966 


1 


-220,943 


2 


37,362 


3 


-6,821 


4 


749 
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Figure 1 1 shows the resulting magnetic field at a series of radii (0, 2.5, 5, 7.5, 
10 and 12.5 cm) plotted in the axial direction produced by the distribution of 
magnetic material, according to the present invention, as shown in Figure 9. Figure 
12 shows the deviation of the above magnetic field from the desired, unifomi one, 
in parts per milUon (ppm), over the ROL Figure 13 shows how the field diverges 
rapidly outside the ROI in the axial direction, but remains reasonably homogeneous 
in the radial direction; see also Table 3 below. 



Table 3: Maximum field deviation in the ROI at each radius. 





Radius cm (% fraction of magnet) 


Maximum deviation in ROI 




0.0 (0) 


0.460 ppm 




2.5 (10) 


0.160 ppm 


m 


5.0 (20) 


-0.541 ppm 




7.5 (30) 


-1.133 ppm 




10.0 (40) 


-1.011 ppm 




12.5 (50) 


0.1 64 ppm 


K U S 

?S ;; 
n 'is? 


15.0 (60) 


2.01 1 ppm 




20.0 (70) 


-4.730 ppm 



As demonstrated by the above figures, the resulting field is unifom to a very high 
degree in the central region. 

Figure 11 serves to highUght the fimdamental inapplicabihty of methods that 
attempt to invert the field equations to the design of permanent magnetic stractures. 
Due to the complexity of the field, it is inconceivable that even someone highly 
conversant with the desi^ of magnets could aibitrarily specify a field such as that 
of Figure 11 with the required degree of accuracy to ensure a practically realizable 
distribution that would yield a reasonably good approximation to the desired field in 
the ROI; essentially this would require a knowledge of the result of the process 
before the process had begun. Certainly, this would not yield the uniformity attained 
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by the present technique which only requires the specification of the known, desired 
field within the ROI. 

In accordance with a second embodiment of the present invention, the design 
method is appHed to homogenize an existing field by distributing magnetic material 
on the curved surface of a circular cylinder, of the type shown in Figure 4. For the 
purpose of illustrating this example, the existing field is assumed to be oriented in 
the z-direction and is substantially uniform except for a first order (linear) 
transverse variation (inhomogeneity) in the ROI. It is to be appreciated that the 
design method is by no means limited to compensating for a linear variation, but is 
generally apphcable to shimming for any inhomogeneity. 

As in the previous embodiment, the magnetic surface is arbitrarily assumed to 
be 0.5 m in diameter and 2 m in length, with a cyUndrical ROI of diameter 0.25 m, 
and length 0.4 m. Optimization was performed on a representative curved 
cylindrical surface of diameter 0.2 m and length 0.4m to produce a transverse linear 
gradient field of 0.1 milliTesla per metre. Again, table 5 below and Figure 17 
demonstrate that this is indeed representative of the ROI. 

For simplicity it will again be assumed that the thickness of the magnetic shell 
is small compared to the diameter of the magnet, for the present purposes 
infinitesimal. In practice, for the creation of shim structures this approximation will 
hold and more elaborate calculations will generally not be necessary. 

From symmetry considerations, an azimuthal variation of cos 
corresponding to a linear transverse variation, in the field can only be produced by 
a first order variation (m=l) in the magnetic material. This again reduces the 
problem to a one-dimensional one of determining just the axial variation. 

The FT of the field can be determined directly fi-om equation 5a. Optimization 
of the shim is then conveniently solved by least squares minimization of the first 
few terms (Cn) over the ROI in real, as against kernel, space. 

In this embodiment coefficients 1, 2, 3 and 4 were necessary as can be seen 
firom table 4. From inspection, it was noted that the zeroth coefficient contributed 
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negligibly to the central field, so this too could be omitted. The amplitude of the 
required axial distribution is shown in Figure 15. The distribution varies azimuthally 
as cos(e) (with amplitude = 1.28 x 10* amp/m) on the surface of the cylindrical 
structure to provide a magnetic field that varies optimally, hnearly, transversely 
5 (Bz(x) = X. 1 00 microTesla/m) over the ROI. 

Figure 16 shows how a shim magnetic stracture designed in accordance with 
the invention and using diamagnetic and paramagnetic materials as illustrated in 
Figure 14 may be implemented in practice, to generate the distribution shown in 
Figure 15. 

0 Figure 17 shows the resulting gradient magnetic field at a series of radii (0, 

2.5, 5, 7.5, 10 and 12.5 cm, 15 cm) , plotted in the axial direction, produced by the 
distribution of magnetic material shown in Figure 15. Note that the Unearity remains 
excellent far outside the ROI. 



Table 4: The harmonic amplitudes for a linear transverse field. 



Harmonic No. 


Amplitude (amps/m) 


1 


247 


2 


-59 


3 


15 


4 


-2 



The resulting magnetic field is linear to a very high degree in the central 
region, see table 5 and figure 17. Certainly, the degree of linearity is well within the 
limitations of first order approximation adopted and of the errors associated with 
manufacturing. Consequently, such a shim designed in accordance with the present 
invention will correct the linear inhomogeneity to a very high order of accuracy. As 
mentioned above, the design method of the present invention is not limited to 
correction of linear variations, and a similar to degree of correction can be 
anticipated generally, for any conceivable inhomogeneity. 
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Table 5: Center field in the ROI at each radius for a linear transverse shim. 



IVitUIUa vlll ^ /O IF AvlIUil UI UId.gUCl^ 


¥t7 tf^micronTMla^ 


0.0 ( 0) 


0.000000 


2.5 (10) 


2.499998 


5.0 (20) 


4.999995 


7.5 (30) 


7.499989 


10.0 (40) 


9.999974 


12.5 (50) 


12.499929 



It is to be appreciated that, generally, the lowest spatial frequency harmonics 
5 are chosen, but this will not always be necessary as witnessed by the example 
above. 

The number of terms over which the optimization is performed may be 
reduced by incorporation of an additional, hypothetical, arbitrary distribution of 
magnetic material over the surface of the structure into the set of spatial frequency 
10 harmonics. By "arbitrary" it is to be understood that a suitable distribution is 
determined by appropriate consideration of the desired field. For example, to 
achieve a field that varies linearly along a flat surface of magnetic material, an 
appropriate distribution might be one whose magnetization also varies linearly 
along that surface. 

15 The method of the invention can be applied to the design of shims in several 

ways. One of these is to fit the dipole harmonics directly to the entire field 
inhomogeneities, as illustrated above. Another is to fit the dipole harmonics to 
individual spherical harmonic terms of the field. The rationale for the latter is that it 
is common practice to map the field in these linearly independent spherical 

20 harmonic terms. The required shim is then the linear combinations of the shims for 
the individual spherical harmonic temis weighted appropriately. This second 
method is advantageous when shimming a series of magnets of the same physical 
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dimensions, as the shim dimensions remain constant and consequently the spherical 
haraionic shims forai a standard set. The resulting shim from the example above 
effectively also represents one component of such a set. 

The required distribution of magnetization can be achieved in several ways. In 
one preferred embodiment the thickness of a single magnetic material can be varied 
in proportion to the required distribution. This is perfectly adequate for the 
production of uniform fields such as described in the first embodiment provided the 
thickness of the material is small in comparison with the sizes of the structure and 
the ROL Where this is not the case the process must be refined by shimming along 
the lines outlined for the second embodiment. 

For tiie purpose of shimming the use of a single unidirectional material cannot 
produce a magnetization that varies sinusoidally, but if an additional uniform field 
can be tolerated simply varying the thickness of the material sinusoidally in the 
manner shown in Figure 18 produces acceptable permissible harmonic 
magnetization components of the form A-cos(m0) and A-sin(m0), in the example 
illustrated, 2-sin(3e). 

In another preferred embodiment distributing materials of differing 
magnetization in the vicinity of the ROI can be used to achieve the desired local 
magnetization. By using differing thicknesses of materials that have opposite 
susceptibilities harmonic components of the form cos(m9) and sin(m9) are possible. 
For example, for a cylindrical surface consisting of a semicircular (0 < 9 < 7i) 
distribution of paramagnetic aluminium (xm = 2.3x10'^), whose thickness varies as 
0.435 A |sin(9)| and an opposing semicircular (-7t<0<O) distribution of 
diamagnetic copper (xm = -10x10'^) whose thickness varies as A |sin(0)|, produces a 
magnetization that varies as A sin(9) to first order ower the whole surface 
(-7C < 9 < 7i), in a similar manner to that shown in Figure 14. Altematively, the 
magnetization distribution shown in Figure 15 could be achieved by machining the 
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above sinusoidal contour on cylinders of paramagnetic and diamagnetic material in 
this manner. 

A material of any required magnetism can be constructed by combining 
paramagnetic and diamagnetic materials in the required proportions. For example, a 
5 material that has effectively zero magnetization can be formed from a combination 
of 303% alimiiniimi and 69.7% copper. This could be a great value in 
manufacturing objects for use in magnetic fields that would not affect the field, or 
be affected by the field. In practice, the proportions may have to be altered slightly 
to take into account the higher order effects when such materials are combined. 
10 For either creating a uniform field or altering an existing one, the 

thickness of the magnet will be a significant fraction of the diameter and the 
consequent placement of material will affect the actual attainable field. To take 
account of this, the method can be applied iteratively as follows, by: 

-calculating the distribution of magnetic material to produce an optimum 
15 approximation to a desired magnetic field in the ROI, in accordance with the 
method described herein 

-determining the total resultant field in the ROI taking into account the 
thicknesses and placements of magnetic materials involved 

-calculating a corrected magnetic material distribution representative of the 
20 difference between said total resultant field and said desired field 

-repeating the latter two steps until the field attained approximates the 
desired one in the ROI to a sufficient degree. 

It is to be appreciated that the homogeneity of the field generated by a 
magnetic structure designed in accordance with the present invention can be further 
25 enhanced by additional active shimming. Resistive electromagnetic shimming may 
subsequently be employed to further improve the homogeneity during the use of the 
magnet. This is particularly relevant at high field strengths with gradient echo 
techniques where the natural susceptibility of the sample distorts the field, in other 
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words, where resistive electromagnetic shimming can be employed to improve the 
homogeneity for each sample. 

Whilst the above embodiments have been described in the context of magnets 
of cylindrical form, it is appreciated that the method of the invention is equally 
5 applicable to the design of magnets of oflier geometric forms, one example of which 
is the parallel flat plate arrangement shown in Figure 1 . 
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CLAIMS. 

L A method of designing a peimanent magnetic structure for generating a 
desired magnetic field in a region of interest (ROI) comprising: 
5 -setting the dimensions of the permanent magnetic structure 

-setting a hypothetical distribution of magnetic material over the surface of 
said structure for each member of a finite set of spatial frequency harmonics of unit 
magnetisation amplitude 

-calculating the magnetic field in the ROI firom each member of said set 
10 -calculating amplitudes of said harmonics by the method of least squares in 

order to achieve the optimum approximation to the desired field in the ROI 

-scaling the members of said set according to said optimised harmonic 
amplitudes 

-establishing, by summing said scaled members, the required distribution of 
15 magnetic material in order to generate the desired magnetic field. 

2. A method according to claim 1 wherein the desired field represents a 
modification of an existing magnetic field. 

3. A method according to claim 2 wherein the desired field represents the 
deviation &om a required field generated by an existing magnetic structure. 

20 4. A method according to claim 3 wherein the desired field represents one of 

the harmonics of the magnetic field over the ROL 

5. A method according to claim 4 wherein each said harmonic results in a 
shim, to create a set of shims representative of the weighted simi of the individual 
harmonics. 

25 6. A method according to any of claims 1 to 5 wherein the desired field is 

produced by a distribution of magnetic material over one or more flat surfaces. 

7, A method according to any of claims 1 to 5 wherein the desired field is 
produced by a distribution of magnetic material over an elliptical cylindrical shell. 
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8. A method according to any of claims 1 to 5 wherein the desired field is 
produced by a distribution of magnetic material over the surface of a circular 
cylindrical shell. 

9. A method according to claim 7 or claim 8 wherein the thickness of 
5 magnetic material is varied in an azimuthal direction in accordance with the desired 

field. 

10. A method according to claim 1 wherein the desired field is produced by a 
combination of paramagnetic and diamagnetic materials. 

1 1. A method of designing a permanent magnetic structure for generating a 
10 desired magnetic field in a region of interest (ROI) on the basis of the minimization of 

the sum of the least squares of the deviations over the ROI. 

12. A method according to any of claims 1-11 wherein an additional, 
hypothetical, arbitrary distribution of magnetic material over said surface is 
incorporated into said set 

15 13. A method according to claims 2 or 3 comprising: 

-calculating the distribution of magnetic material to produce an optimimi 
approximation to a desired magnetic field in the ROI, in accordance with the method 
described herein 

-determining the total resultant field in the ROI taking into account the 
20 thicknesses and placements of magnetic materials involved 

-calculating a corrected magnetic material distribution representative of the 
difference between said total resultant field and said desired field 

-repeating the latter two steps until the field attained approximates the desired 
one in the ROI to a sufficient degree. 
25 14. A method according to any of claims 1 to 5 where the magnetic fields of 

said set of harmonics are evaluated by an actual step involving numerical integration 
in real space. 
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15. A method according to any of claims 1 to 5 where the Fourier 
transformers of the magnetic fields of said set of harmonics are evaluated directly in 
Fourier Kernel space. 

16. A permanent magnetic structure designed in accordance with the method 
5 of any of claims 1 to 15. 
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Fig.6. 
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Fig.8. 
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Fig. 10a. 



Fig. 10b. 
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Fig.11. 

Permanent Magnet Uniform Field 
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Fig. 12. 

Permanent Magnet Uniform Field 
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Fig. 13. 

Permanent Magnet Uniform Field 
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Fig. 17. 

Passive linear transverse gradient field 
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